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Abstract —An alternative parametric description for discrete 
random variables, called muculants, is proposed. Contrary to 
cumulants, muculants are based on the Fourier series expansion 
rather than on the Taylor series expansion of the logarithm of 
the characteristic function. Utilizing results from cepstral theory, 
elementary properties of muculants are derived. A connection be¬ 
tween muculants and cumulants is developed, and the muculants 
of selected discrete random variables are presented. Specifically, 
it is shown that the Poisson distribution is the only discrete 
distribution where only the first two muculants are non-zero, thus 
being the muculant-counterpart of the simple cumulant structure 
of Gaussian distributions. 

Index Terms —Higher-order statistics, cepstrum, cumulants, 
discrete random variables, Fourier series expansion, Edgeworth 
expansion, homomorphic systems 


I. Introduction 

Cumulants are parametric descriptors of random variables 
(RVs) - no matter whether continuous or discrete - and are 
commonly used to analyze non-Gaussian processes or the 
effects of nonlinear systems (TJ. Unlike moments, cumulants 
are orthogonal descriptors and satisfy a homomorphism prop¬ 
erty 12, 0: The cumulants of a sum of two independent RVs 
are the sums of the respective cumulants of the individual RVs. 

As attractive this additivity of cumulants is, there are several 
drawbacks in approximating an RV just in terms of these 
descriptors: First of all, since the cumulants are the Taylor 
series coefficients of the logarithm of the distribution’s char¬ 
acteristic function 0, they constitute a complete description of 
the distribution only if the characteristic function is infinitely 
differentiable. Secondly, even if all cumulants exist, there are 
no general results regarding the behavior of the sequence 
of cumulants; in fact, the sequence may even diverge: For 
example, the cumulants of the continuous uniform distribution 
are the Bernoulli numbers, a sequence which is known to 
diverge. In such a case, also approximating the distribution 
or characteristic function of an RV based on finitely many 
cumulants is problematic, since the divergence of the cumulant 
series prevents bounding the approximation error. Addition¬ 
ally, according to the Marcinkiewicz theorem, the Gaussian 
distribution is the only distribution which can be characterized 
with finitely many (two) cumulants: Every approximation of a 
non-Gaussian distribution with other than two cumulants will 
fail to correspond to a distribution at all. It also follows that 
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cumulants are not well suited for hypothesis testing, except in 
the important case of testing the Gaussianity of a distribution. 

In this paper, we address these problems by replacing the 
Taylor series expansion of the logarithm of the characteristic 
function by a Fourier transform. This approach guarantees that, 
provided the Fourier transform exists, the approximation error 
for a truncated integration range can be bounded, and that this 
bound converges to zero for increasing integration ranges. This 
benefit is bought at the cost of a functional (i.e., uncountable) 
description rather than a parametric (i.e., finite or countable) 
one. 

The situation greatly simplifies when we place the focus 
on discrete distributions supported on the set of integers. For 
these, the characteristic function is periodic, and the Fourier 
transform can be replaced by a Fourier series expansion. 
This approach is not only the most natural one (periodic 
functions lend themselves better to Fourier analysis than to 
a Taylor series expansion), but the Fourier series coefficients 
- henceforth called muculants - also provide a parametric 
description of the distribution of the RV. Hence, in this work, 
we take a first step toward an alternative description for RVs by 
focusing on discrete distributions only. Specifically, we show 
that the Poisson distribution is the only distribution for which 
only the first two muculants are not zero. In this sense, the 
Poisson distribution is the Fourier-analysis counterpart of the 
Gaussian distribution, which is characterized by having only 
two non-zero Taylor series coefficients. This result opens a 
new way to test whether a distribution is Poisson, just as 
cumulants are analyzed to test Gaussianity. 

The sequence of operations: Fourier transform - (complex) 
logarithm - inverse Fourier transform, that needs to be taken 
to obtain the characteristic function from a distribution, and 
to obtain the Fourier series coefficients from the logarithm 
of the characteristic function, is also essential in cepstral 
analysis 0, 0. In this sense, our work follows the spirit 
of Widrow, who also applied methods from signal processing 
to the analysis of RVs 0, 0. Specifically, he argued that the 
nonlinear quantization operation is a linear sampling operation 
of the probability density function of the quantizer input 
signal. This allowed him to propose a quantization theorem 
reminiscent of Shannon’s sampling theorem, stating conditions 
under which the distribution of the quantizer output can be 
used to reconstruct the distribution of its input. We show 
that taking this signal processing-inspired approach one step 
further can provide new insights into the theory of discrete 
distributions. Whether cepstral analysis can also be applied to 
analyze operations on RVs (such as quantization), is within 
the scope of future investigation. 
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A few words regarding the naming are now in order: 
The name cepstmm is derived from a reversal of the first 
syllable of the word spectrum : While the cepstrum - being 
the inverse Fourier transform of a nonlinear operation of 
a Fourier transform - is actually situated in the original 
(e.g., time) domain, Bogart et al. introduced new terminology 
to emphasize that this sequence of operations can provide 
fundamentally different insights than would be possible by 
thinking in the original domain only. For example, they got 
frequency variables in the first, and quefrency (rather than 
time) variables in the second Fourier transform, and called 
a convolution (filtering in time domain) liftering in cepstral 
domain. Following this approach, we decided to call our new 
parametric descriptors muculants, a reversal of the first syllable 
of the word cumulants. 


A. Related Work 

Theory regarding the characteristic function, its behavior 
and some of its applications, can be found in the work 
of Eugene Lukacs a, 0 , ca. For elementary results on 
cumulants the reader is referred to im and m. 

To reconstruct a distribution function in terms of its cu¬ 
mulants, the Edgeworth-expansion is well known fl3l , fl4l . 
Besides that, there exist expansions applying orthogonal poly¬ 
nomials, see EQ and the references therein. Nevertheless, 
such approximations based on cumulants or moments need not 
converge for arbitrary distributions, neither need a truncated 
approximation correspond to a valid distribution. Hence, if 
we want to describe non-Gaussian distributions, it might be 
promising to look for new descriptors. A distribution can be 
fully characterized in terms of its moments if all moments 
exist, respectively if the cumulant series does not diverge. 
Damping the exponential function to guarantee this conver¬ 
gence led to the definition of q-moments and q-cumulants 

in css, ca. 

The Marcinkiewicz theorem 0 states that any function, 
which is the exponential of a polynomial and is positive 
definite on R, has a polynomial which is of at most degree 
two. Since the cumulants are obtained from the derivatives of 
the logarithm of the characteristic function, this means that 
any RV either has infinitely many cumulants or has cumulants 
up to second order only. Clearly, every approximation which 
is based on cumulants can only take finitely many coefficients 
into account. Hence, the reconstruction is either Gaussian or 
does not correspond to a RV at all. 

The cepstrum was introduced and first published in 0. The 
main motivation was to investigate the influence of an echo 
in a signal and how the echo delay can be determined by 
looking at the cepstrum. Around the same time, Oppenheim 
introduced similar ideas in his dissertation M under the more 
general motivation of representing nonlinear systems in a way 
to ensure the principle of superposition. Today the cepstrum is 
widely used in speech processing Q. The whole terminology 
of cepstral analysis is closely related to these applications, 
therefore we shall only refer to the term cepstrum if absolutely 
necessary to avoid unnecessary confusion. 


B. Organization of the Paper 

Based on the definition of the cepstrum and the properties 
of the characteristic function, we present two definitions for 
muculants in Section III-AI and discuss their properties in 
Section III-CI Utilizing results from signal processing, the 
existence of muculants is related to the Paley-Wiener theorem 
in Section IH-B1 while the connection between the two com¬ 
peting definitions is shown to be connected to minimum-phase 
systems in Section IH-DI The connection between muculants 
and cumulants is discussed in Section [Till Finally, Section HVl 
is devoted to deriving the muculants for selected distributions, 
such as Poisson, Bernoulli, geometric and binomial. 

IT Definitions and Properties 
Let X be a real-valued RV, and let T 1: (/i) denote its 
characteristic function 

**(/i) :=E(e^ x ). (1) 

The characteristic function always exists since it is an expec¬ 
tation of a bounded function. There is a direct correspondence 
between <b T (//) and the distribution of X: Two RVs are iden¬ 
tically distributed if and only if their characteristic functions 
are equal. In addition to that, a characteristic function has the 
following properties relevant for our work C3p. 13]: 

1) <!>.,. (//) is uniformly continuous everywhere. 

2) $a:(M)U=0 = l- 

3) $x(m) is bounded. On its entire domain, |$ x (/z)| < 1. 

4) $x(m) is Hermitian. 

Finally, let X and Y be two independent RVs and let Z = 
X + Y be their sum. Then, the characteristic function of Z is 
given as 

<M M ) = $ X (p) ■ $y(p). (2) 

Letting X be a discrete RV taking values from the set of 
integers Z, it can be described by its probability mass function 
(PMF) f x : Z —>• [0,1], where 

<E Z : / x [f] :=P(X = 0- (3) 

In this special case, the characteristic function is the inverse 
Fourier transform of the PMF and periodic, as becomes 
apparent from 

OO 

<m^)= e w 

—OO 

Moreover, letting X and Y be two independent, discrete RVs 
with PMFs f x and f y , respectively, the PMF of their sum, 
Z = X + Y, equals the discrete convolution of the respective 
PMFs, i.e., 

fz = fx * fy (5) 

We call a PMF causal if f x [£\ = 0 for £ < 0. We call a 
PMF minimum-phase if it is causal and its ^-transform 

OO 

^x{z) := E-^]^ (6) 

S=o 

has all its zeros inside the unit circle (possible poles are inside 
the unit circle by construction, since every PMF is absolutely 
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summable). Note that the characteristic function is related to 
this z-transform via 

^ x {p)=^ x {e^). (7) 

A. Types of Muculants 

We define two types of muculants, analyze their properties, 
and develop how they are connected. The definitions mainly 
follow the definitions of the cepstrum given in lf20l Ch. 13], 
with the main difference that the roles of Fourier transform 
and inverse Fourier transform are reversed (the characteristic 
function is the inverse Fourier transform of the PMF, rather 
than the Fourier transform). While some properties of the mu¬ 
culants are directly transferred from cepstral analysis, several 
of the presented results are based on the fact that we operate 
on a probability space. 


Since the natural logarithm is unambiguous by definition, 
the power muculants can be computed by 

c 4 n ] = ^fJ ^\<$> x {p)\ 2 e~ mn dp (11) 

whenever the integral exists. By construction, the power mu¬ 
culants cannot be used to compute the characteristic function 
of a RV in general, since the magnitude function is not 
invertible. Nevertheless, under specific conditions the phase 
information of the characteristic function can be retrieved 
from its (squared) magnitude. We will elaborate on this in 
Section IH-DI 

It is worth mentioning that the Mercator series for complex 

z, 

log(l + z)=V^ -(12) 

z ' m 

m —1 


Definition II. 1 (Complex Muculants). The complex mucu¬ 
lants {c x [n] } n( zz are the coefficients of the Fourier series 
expansion of logT^/i) (if it exists), i.e.. 


log$ x (p) = c x [n\e mn . (8) 

n——oo 

where log denotes the complex logarithm. 

Note that the phase of the complex logarithm - defined 
as the inverse operation of the complex exponential - is not 
uniquely defined since an integer multiple of 27 t can be added 
to the phase term without changing the evaluation. Hence, the 
computation of the complex muculants requires a continuous 
phase | 20l p. 1009]. Assuming that log<l> x (/u) has no zeros on 
the unit circl^l we can assign such a continuous phase and 
unambiguously define the complex logarithm as CD 


log z = \n \z\ + j ■ arg*(z), (9) 


where In is the natural logarithm and where arg*(z) = 
Arg(z) + 2kn for z € C and k £ Z such that phase jumps 
by 27 t are resolved. One can formally replace the z-plane 
with a carefully defined Riemann surface, and thus represent 
log(z) as a single-valued analytic function |]22l . Since the 
computation of such an “unwrapped” (i.e., continuous) phase 
function is essential for the proper estimation of the complex 
cepstrum, a vast body of literature can be found to deal with 
this problem in practical cases |20l . |2P . |f23l , 1241 . 

Letting log(z) be unambiguously defined in the above sense, 
we can estimate the complex muculants through 


AA = ^f log^(/x)e-^" dp (10) 


whenever the integral exists (cf. Section Hl-Bl i. 

In case we are only interested in the magnitude of the char¬ 
acteristic function, we may resort to the power muculants: 


Definition II.2 (Power Muculants). The power muculants 

{c x \ri] }n(? 7 /, are the coefficients of the Fourier series expansion 
of In |<b; r (/i)| 2 (if it exists). 


where \z\ < 1, z ^ —1, can be used to derive a recursive 
relation among the complex muculants, at least under special 
conditions. Such a relation, obviating the need for integration, 
has also been described in ll20l p. 1022] and might be of 
interest for numerically computing the muculants for PMFs 
not given in explicit form. 

In what follows, we will present properties for complex 
muculants only. By using their relation to power muculants 
introduced in Section Hl-DI below, most of these properties can 
be carried over to them as well. 


B. Existence of the Muculants 

In this section, we present a sufficient condition for the ex¬ 
istence of muculants; the discussion of whether the muculants 
can be used to describe the characteristic function is deferred 
to Proposition III. 21 in Section ITl-CI 

The Fourier series expansion of a function exists if the 
function is absolutely integrable l25l p. 402]; hence, the 
complex muculants (see Definition III. II ) exist if 


1 

27T 


I log$ x (p)\dp < oo. 


Using (O and applying the triangle inequality, we get 


(13) 


I In I'M/')!! < |log 0Mm))| < |ln|$z(M)|| + |arg*($ x (/z))|. 

(14) 

The phase of the A> x (//.) is a continuous function (cf. properties 
of the characteristic function) with compact support [— 7r, 7 t], 
thus its integral is finite. Absolute integrability is consequently 
guaranteed if and only if 


1 

27T 


I In \$> x (p)\\dp < oo. 


(15) 


It is apparent that this condition for existence of the com¬ 
plex muculants also applies to the power muculants. Since 
I'Mfi)! < 1 and, thus, log < 0, we can rewrite <m 

as 


1 

27T 


In \<& x (p)\dp > — oo. 


(16) 


'If t-fpi) has zeros on the unit circle, phase jumps by rr are introduced Note that these condition rules OUt the existence of muculants 

and as a consequence one cannot determine an unambiguous phase curve (Til for characteristic functions that vanish on an interval. 
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The condition in ( IT6l ) reminds of the Paley-Wiener condi¬ 
tion |26l p. 423]. The Paley-Wiener condition, loosely speak¬ 
ing and translated to the language of probabilities, states that 
to the given magnitude a unique phase can be associated which 
guarantees that the corresponding PMF is minimum-phase. 
Hence, if we know a priori that the PMF is minimum-phase, 
then the characteristic function can be described completely 
via the power muculants (see Section IH-DI below). Interest¬ 
ingly, causality alone is not enough: Not every causal PMF 
is minimum-phase, and hence, not for every causal PMF the 
power cumulants are a complete description l27l p. 215]. 

C. Properties of the Muculants 

This section will list some properties of the muculants 
resulting from the fact that we investigate PMFs. Throughout 
this section we assume that the muculants exists. 

Proposition II. 1 (Properties of the Muculants). The complex 
muculants satisfy the following conditions: 

1) c x [n] e R for all n; specifically, if f x (f) = f x (—Q, 
then Cx[n] = c x [—n] 

2) c x [ 0] < 0 

3) lim c x [n\ = 0 

n—f^zoo 

Moreover, if X and Y are independent RVs, then the mucu¬ 
lants of their sum, Z = X + Y, satisfy 

c z [n] = c x [n\ + c y [n]. (17) 

Proof: The first property follows from the fact that the 
characteristic function ( l> : , (p) is Hermitian, and from the fact 
that the Fourier transform of a Hermitian function is real¬ 
valued |20l p. 83]. Specifically, if the PMF is an even function, 
then the characteristic function is real and the muculants have 
even symmetry. 

For the second property, note that 

4[0] = ±J log $ x (p)dp (18) 

= j In I ^ x (p)\dp (19) 

( b) 

< 0 ( 20 ) 

where (a) follows from (0 and the fact that the phase of the 
characteristic function has odd symmetry]! (J>) then follows 
from |<b a; (/f)| < 1. 

The third property is a direct consequence of the Lebesgue- 
Riemann theorem and the fact that log$ x (^) is absolutely 
integrable ll28l p. 94]. 

Finally, if X and Y are independent RVs, then 4> x+y (//) = 
$ x (/r) • Qy(p). The complex logarithm with its unambiguous 
phase defined in 0 satisfies]! 

log(w^) = log(wj) + log(^) (21) 

= In |iu| + In \z\ + j ■ arg *(w + z ) (22) 

= ln|u;| + j • arg *(w) + ln|z| + j • arg* (2), (23) 


for w, z e C. The desired result follows. ■ 

That muculants are additive descriptors, i.e., that the mu¬ 
culants of a sum of independent RVs are the sum of the 
respective muculants, was one of the main motivations for 
introducing homomorphic transformations in first hand: The 
benefits of cumulants are retained while eliminating some of 
their drawbacks particularly problematic with discrete RVs. 

Note that, for the power muculants, |<f> x (/r)| 2 is the char¬ 
acteristic function of the difference between two independent, 
identically distributed (i.i.d.) RVs with characteristic functions 
& x (p). It follows that the PMF corresponding to |$ ;c (/x)| 2 
equals the autocorrelation function of the PMF f x 126] p. 182]. 
Since the autocorrelation function of a real-valued function 
has even symmetry, from the first property of Proposition III. II 
follows that the power muculants are always even symmetric. 

It remains to show in which case the complex muculants 
are sufficient to describe the characteristic function of a RV. 

Proposition II.2 (Convergence of the Muculants). Let the RV 
X have finite expectation. Then, the Fourier series expansion 
of\og<V x (p) - if it exists - converges pointwise, i.e., 

OO 

log$x(/T)= £ x\n\ emn , (24) 


Proof: The Fourier transform for the 27T-periodic function 
log$3;(/x) converges pointwise if log<F x (/r) is piecewise con¬ 
tinuous and piecewise differentiable ll28l p. 105]. Piecewise 
continuity of log<l> x (/r) is guaranteed by uniform continuity 
of '[> , (//) and piecewise continuity of log. For differentiability, 
note that 


d_ 

dp 


log <V x (p) 


$x(pY 

*&x(p) 


(25) 


where (p) can vanish on an at most countable set of points 
(otherwise the muculants would not exist, cf. Section III-BI) . 
Finite expectation guarantees existence of the first derivative 
of $ x {p) \m P- 4^1’ which completes the proof. ■ 

Note that we cannot drop the assumption of finite expecta¬ 
tion, because one can construct characteristic functions which 
are nowhere differentiable. For example, the Weierstrass func¬ 
tion <& x (p) = £^=1 (5) cos(3 n p) is nowhere differentiable, 
but continuous everywhere and a valid characteristic function, 
see Jl9] p. 47]. 

We can strengthen the third property of Proposition III. II if 
the muculant expansion converges: 


Proposition II.3. Let the muculant expansion converge as in 
Proposition I//.2I Then, c x [n] = 0(1/n). 


Proof: The conditions for pointwise convergence imply 
bounded variation l29l p.70]; the result about the order of 
convergence follows from (30). ■ 


Proposition II.4 (Sum of the Coefficients). If the RV X has 
finite expectation, then 


2 Even if we would allow phase jumps by n, we require that the character¬ 
istic function is Hermitian and thus the phase term vanishes. 

3 Note that this property does not hold for the principal value of the complex 
logarithm mo P . ioo9]. 


^2c x [n\ = 0 

n 

provided the muculants exists. 


(26) 
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Proof: Finite expectation of X guarantees that 

OO 

log$a;(/Lt)= ^2 c x [n]e mn (27) 

n =—oo 

pointwise. Evaluating the left-hand side at p, = 0 yields zero, 
since < f , x (At)|/ J =o = 1. For p = 0 the right-hand side equals 
the sum of the muculants, which completes the proof. ■ 
The fact that muculants must sum to zero makes truncating 
the Fourier series to approximate a distribution as problem¬ 
atic as truncating the cumulant expansion: The truncated 
series need not correspond to a valid PMF. The fundamen¬ 
tal difference, however, is that by Parseval’s theorem and 
Proposition III. 31 the squared error between the true and the 
approximated characteristic function does not diverge but stays 
bounded. Hence, muculants may be superior to cumulants 
when used in functionals of distributions, such as entropy 
or informational divergence. It would be interesting to in¬ 
vestigate, e.g., if an approximation of negentropy based on 
muculants is possible and has better convergence properties 
than one based on cumulants EE 

D. Connection Between the Complex Muculants and the 
Power Muculants 

The connection between the different types of muculants is 
based on ©; this equality guarantees that log z* = (logz)* 
and that with log \z\ 2 = log z + log z* one has 


magnitudes of <f> x (/r) and $x,min(d) are equal, the minimum- 
phase assumption implies that the complex muculants of the 
minimum-phase component are given by the power muculants 
of the original distribution via d29l . i.e., 

(c x [0 ]/2, if n = 0 

— \ if Tl 0 • (31) 

[o, else 

Since the convolution in equation (l30l > resolves to an addi¬ 
tion of the muculants, the allpass contribution is given by 

Cy — C x C x ,min- ( 32 ) 

The corresponding PMF can be obtained by applying the ex¬ 
ponential function to ([8j and evaluating the Fourier transform 
subsequently. Such a decomposition clearly works for PMFs 
of RVs which are bounded from below, i.e., if X > M. In 
this case the allpass component is a degenerate distribution 
with PMF /*,[£] = <$[£ — M], The minimum-phase component 
is thus an RV which is non-negative, and which has its first 
mass point in the origin. 

That such a decomposition can also fail to deliver two valid 
PMFs is shown in the following example: 

Example II.l. Fet f x be the PMF of a left-hand sided 
geometric distribution with p = 0.2, such that: 


c x [n] = c x [ n ] + c x [~n]. (28) 

according to the Fourier transform properties l20l p. 83]. 
We can therefore directly obtain the power muculants by 
adding pairs of complex muculants. Such a correspondence 
is obviously not reversible in general. However, considering 
the concept of minimum-phase sequences |20j p. 1013], we 
obtain conditions under which it is possible to derive complex 
muculants from power muculants. A key result in this regard 
is that the complex cepstrum of a minimum-phase sequence is 
causal l20l p. 1013]. As such, if f x is minimum-phase, then 
c x [n\ = 0 for n < 0, making ( 128l > reversible: The complex 
muculants can be computed from the power muculants by 


/*[£] 


p(l -p) if £ < 0 

0 , if £ > 0. 


(33) 


Applying above decomposition, it can be shown that f x , m in 
is a right-hand sided geometric distribution, i.e.. 


fx^rain [£] 


0 , if £ < 0 

if£>0 


(34) 


This is indeed a valid PMF. However, it can be seen that f y 
violates the constraint for a PMF to be non-negative. Actually 
the allpass component, given by 


I c x [0]/2, if n = 0 

c x [n] = < c x [ri], if n > 0 . (29) 

[o, else 

Hence, the minimum-phase property implies a relationship 
between phase and magnitude of the Fourier transform, known 
as the Hilbert transform. Having such a relationship, it is 
sufficient to evaluate the power muculants. 

The question arises if we can decompose every given PMF 
f x into two parts, one of which has the minimum-phase 
property: 

fx[t] = (U, m in*f V )[Z\ (30) 

Such a decomposition is always possible, as long as the com¬ 
plex muculants exists (see Chapter 5.6.1 on minimum-phase 
and allpass decomposition in lt20l ). but it is not immediately 
clear whether both f x ,min and fy are valid PMFs. Since the 


(p{2-p)(l - k) t, if £ < 0 
/»[£] = <-(1-p), if £ = 1, (35) 

[o, else 

is negative for £ = 1. 

Interestingly, however, f y [£] = I is fulfilled; this is true 
in all cases in which the minimum-phase component is a valid 
PMF, since by (f30t 

1= E^1 = EE ^ C] fy [C] (36) 

« « c 

= E fx,min[f, — C]/y[C] = fy[C\- (37) 
C 5 C 

Whether or not the minimum-phase component is always a 
valid PMF is subject of future investigation. 
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III. Linking Cumulants and Muculants 

Using the 2-transform instead of the Fourier transform, a 
close connection between the cepstrum and the cumulants is 
pointed out in J32). In 11331 , the moments of the cepstrum are 
calculated and it is shown how they correspond to the moments 
of the original sequence, in order to circumvent the explicit 
calculation of the cepstrum. 

Motivated by the works mentioned above, we will use basic 
Fourier transform properties to derive an explicit connection 
between the muculants and the cumulants. We start with the 
inverse Fourier transform as relation between the characteristic 
function and the muculants 

OO 

log<I> x (/r) = ^ c x [n]e mn . (38) 

71 =— OO 

Assume that log <I>, (/./.) is periodic and continuous (not only 
piecewise continuous), that its first (fc — 1) derivatives are 
continuous, and that its fc-th and (fc + l)-th derivatives are 
piecewise continuous. Then, ll28l Theorem 3.22] ensures that 
one can represent the fc-th derivative of log <3? x (/x) by the 
following Fourier series: 

7/c 00 

— \og<f> x (p) = J2 (jn) k c x [n\e j ^ n (39) 

' TI— — 00 

Evaluating the left-hand side at p = 0 yields j k rtk, where Kk 
is the fc-th cumulant; evaluating also the right-hand side at 
zero one thus obtains 

OO 

Kk = ^2 n k c x [n\. (40) 

71 =— OO 

The sum on right-hand side of d40l) resembles the fc-th 
moment of a discrete distribution. Abusing terminology by 
ignoring the fact that the muculants need not be non-negative, 
we claim that the fc-th cumulant (if it exists), equals the 
k-th non-central moment of the complex muculants. Note 
further that every single cumulant contains information of all 
muculants. It shall be noted that there is a similar expression 
called Good’s formula 041 . which expresses the cumulants in 
terms of the moments of a random vector. Moreover, there 
is some interest in estimating the cumulants as combination 
of moments m. The question arises if the relation between 
cumulants and muculants provides insight into this problem. 

IV. Muculants of Selected RVs 

In this section we will present the complex muculants for a 
selection of discrete distributions. While the majority of these 
results has been obtained by applying the Mercator series to 
log {p) and by equating the coefficients subsequently, the 
Poisson distribution allows a simpler treatment, as outlined 
below. The complex muculants for all considered distributions 
are summarized in Table Q] 

A. Poisson Distribution 

Letting A G R+, the characteristic function of the Poisson 
distribution is given as 

<S> x (p) = e x ^~ 1 \ (41) 


All its cumulants exist and equal A. A fortiori, X has finite 
expectation and we can write with Proposition III. 21 

OO 

log(e A ^ eJM_1 ^) = Ae JM — A = c x [n] ■ e^ n . (42) 

71 = — OO 

By equating the coefficients we obtain 

{ —A, if n = 0 

A, if n = 1 . (43) 

0, else 

Hence, all but two coefficients are equal to zero. These two 
coefficients directly correspond to the parameter of the Poisson 
distribution. Considering the fact that the sum over all cepstral 
coefficients is always zero and that c x [0] is always non-positive 
(cf. Proposition lH.il ), we can state that the Poisson distribution 
is the only distribution with only the two first muculants non¬ 
vanishing. The Poisson distribution is thus the counterpart of 
the Gaussian distribution, which is unique by being the only 
distribution with just the first two cumulants non-vanishing. 
However, while there are no distributions with finitely many 
(but more than two) non-vanishing cumulants, distributions 
with, e.g., three non-vanishing complex muculants exist. 

B. Degenerate Distribution 

The degenerate distribution is used to characterize constant 
RVs, i.e., those for which X = M holds almost surely. The 
PMF for such a distribution is given by the Kronecker delta, 
i.e., fx[£\ = <5[£ — M]. The characteristic function is given 
by <J> x (/i) = consequently, for the power muculants we 

have c x [n\ = 0. While c x [0] = 0, for n ^ 0 the complex 
muculants are given as 


r , 1 r 


c x[n] = — / jpMe m dp 

(44) 

= — (-l) n+1 . 
n 

(45) 


Note that since in this case log<l> x (gt) = jpM is not periodic 
and continuous (but only piecewise continuous), one cannot 
expect (l40b to hold: And indeed, while K3 = 0, the third non¬ 
central moment of the complex muculants diverges. 

C. Geometric Distribution 

Let X be a geometric distribution with PMF f x [£] = p( 1 — 
p)£ for £ > 0, else f x [f\ = 0. The characteristic function is 
given by <f> x (/x) = i-(i - P y e tr • 

In the integral 

6 - |n| = h /j oge ~" nd "- <46) 

one can split the logarithm and integrate over the numerator 
and denominator separately. For n = 0, the denominator in¬ 
tegrates to zero and, consequently, c x [0] = log p remains. For 
n/0, the numerator vanishes and, with 2 = 1 — (1 — pje - - 774 , 
we can expand log (l — (1 — p) ■ e- 7 ^) with the help of the 
Mercator series given in (IT2l > as 

log (1 - (1 - p) ■ e™) = V (47) 

71=1 
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TABLE I 

Complex muculants c x [ n ] of selected distributions. For comparison, we also present the cumulants n k and the characteristic 

FUNCTIONS 


Distribution 

MZ\ 

'MaO 

Kb 

c x [n\ 

Poisson 

\ 

e- A 4f, £>0 

0, else 

e A(^> —1) 

<< 

II 


—A, if n = 0 

A, if n = 1 

0, else 

Degenerate 

i[f - M 1 


j M, if k = 1 

Kk ~ | 0, else 


0, if n = 0 

else 

Bernoulli 
(p < 0.5) 

■ 

'l-p, if £ = 0 

p, if £ = 1 

0, else 

1 — p + pe J/i 

k„ + i = p(i -p)^r 

1 

Tog(l — p), if n = 0 

(_ T +1 (a)"’ *"><> 

,0, else 

Bernoulli 
(p > 0.5) 


'l-p, if £ = 0 

P, if « = 1 

0, else 

1 — p + pe ik- 

K n+ 1 = P(1 -P)^f 

< 

logp, if n = 0 

^rMvr)- *"<° 

if n > 0 

Geometric 


(1 -p) ( p, £>o 

0, else 

V 

l-(l-p)e-^ 

Kn + l = P(l + P)pfjf 

P = (1 ~P)/P 


logp, if n = 0 

(1-p)n , if n > 0 

0, else 

Negative 

Binomial 

1 

’( e+ e _1 )(i -pYp*, «>o 

0, else 

( i-p V 

r times the cumulant 
of the geometric distribution 


r log(l — p), if n = 0 
rpU if n > 0 

V 1 —peW ) 

0, else 

Binomial 

1 

f(^(i-pf-0, €>0 

0, else 

(l — p + pe-?^) N 

N times the cumulant 
of the Bernoulli distribution 


'lVlog(l— p), if n = 0 

if n. > 0 

,0, else 


The muculants for n > 0 are then obtained via equating the 
coefficients as 

Cl - v) n 

Cx[n\\ n > o =-. (48) 

n 

Interestingly, if the characteristic function (except for a 
constant) can be expanded with the help of the Mercator series, 
it turns out that c x [n] = 0 for n < 0. This corresponds to a 
PMF with minimum phase (cf. Section Hl-Db . 

That this is no coincidence is outlined in the following 
argument. If the PMF is minimum-phase, all poles p k and 
zeros 0^ of the associated 2-transform 


^x(z) 


nr=i(i -Qfc*- 1 ) 

nr=i (i-p^- 1 )’ 


(49) 


lie inside the unit circle, i.e., both |oj.| < 1 and \p k \ < 1 for all 
k. Exploiting & x (p) = (e JAt ) and the logarithmic identities 


logi> x (/i) = ^2 (l - o k z ^-^(l -p k z x ) , (50) 

k= 1 k-1 


one can consequently use the Mercator series to derive the 
muculants. 


V. Discussion and Future Work 


That the Poisson distribution has only two non-vanishing 
muculants makes the presented theory attractive for hypothesis 
testing. Future work shall develop a test whether a distribu¬ 
tion is Poisson or not, which requires numerical methods to 
estimate muculants from data. The recursive relation presented 
in lf20l p. 1022] is of interest in this regard and might facilitate 
the derivation of estimator variances, confidence bounds, and 
optimality statements. 

Closely related is the quest for possible applications of 
muculants. We believe that due to their desirable properties, 
they might yield, e.g., expressions for (neg-)entropy or infor¬ 
mational divergence that converge faster than those based on 
cumulants ED. 

We presented a condition for existence of muculants in 
Section III-B1 however, it is not clear whether there exist 
distributions violating this condition. The search for a con¬ 
crete example, or, preferably, for a more general statement 
about existence is within the scope of future work. Similar 
investigations shall be devoted to the convergence of the 
muculant series (cf. Proposition III. 2I >. which for the moment 
is guaranteed only for distributions with finite first moments. 


While muculants behave better than cumulants when trun¬ 
cating their series in the sense that the resulting error is 
bounded and decreases with increasing number of summands, 
the question under which circumstances a truncated series rep¬ 
resents a valid PMF remains open. Although distributions with 
finitely many (but more than two) non-zero muculants may 
exist, the relation between these distributions corresponding 
(possibly) to a truncated series and the distribution correspond¬ 
ing to the complete series deserves closer investigation. 


Finally, the theory of muculants shall be extended to 
continuous RVs, even though this might require a muculant 
function rather than a muculant sequence. In this context, the 
connection between the muculant function and/or cumulants 
of a continuous RV and the muculants of a discrete RV 
obtained by uniformly quantizing the former might be of great 
interest. A fundamental step towards these results lies in the 
fact that quantization periodically repeats the characteristic 
function 0. 
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